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Abstract 

Large  voltage  step  chronoamperometry  is  shown  to  be  a  time-efficient  means  to  examine  solid  electrode  materials  compared  with  conventional 
electrochemical  methods  such  as  linear  sweep  voltammetry  (LSV)  and  step  potential  electrochemical  spectroscopy  (SPECS),  all  the  while  providing 
comparable  information  concerning  the  rate  capability  of  a  material  and  its  capacity.  The  applicability  of  the  technique  is  demonstrated  through  a 
study  of  the  alkaline  y-Mn02  electrode.  By  sampling  the  current  (and  hence  the  charge)  at  various  times  after  the  chronoamperometric  voltage  step, 
the  compatibility  between  chronoamperometry  and  LSV  is  disclosed.  Furthermore,  modelling  of  the  chronoamperometric  data  using  two  curves 
based  on  a  spherical  diffusion  model  representing  fast  and  slow  discharge  processes  are  found  to  be  statistically  suitable.  From  this  modelling, 
values  of  AJD  (where  A  is  the  electrochemically  active  surface  area  and  D  is  the  diffusion  coefficient)  for  the  two  processes  are  3.89  x  10-4 
and  0.70  x  10-4  cm3  s-1/2  g_1,  respectively,  both  of  which  are  comparable  with  AJD  data  extracted  from  a  SPECS  experiment  on  an  identical 
electrode. 
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1.  Introduction 

In  modern  society  there  is  an  increasing  trend  towards 
electronic  portability  (e.g.,  computers,  compact  disc  players, 
flashlights),  and  as  such  an  increasing  demand  is  being  placed 
on  the  batteries  that  are  used  to  power  these  devices.  The  aque¬ 
ous  alkaline  Zn|MnC>2  cell  is  a  popular  choice  as  a  portable 
power  supply  due  to  the  good  electrochemical  properties  of  the 
manganese  dioxide  cathode.  Despite  this  popularity  and  per¬ 
formance  fundamental  understanding  of  the  behaviour  of  the 
manganese  dioxide  is  very  limited,  due  mainly  to  the  com¬ 
plex  interplay  between  material  structure  and  electrochemical 
performance  [1-3]. 

When  characterizing  the  electrochemical  behaviour  of  a  com¬ 
mercial  solid  electrode  system  it  is  desirable  to  gain  as  much 
information  as  possible  with  a  minimum  number  of  experiments 
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so  as  to  save  sample  quantity  and  be  more  efficient  time-wise. 
There  are  numerous  techniques  available  to  the  electrochemist 
to  characterize  material  performance  including:  (i)  potentiody- 
namic  methods  such  as  linear  sweep  and  cyclic  voltammetry 
(LSV  and  CV,  respectively);  (ii)  potentiostatic  methods  such 
as  chronoamperometry  and  the  more  elaborate  method  of  step 
potential  electrochemical  spectroscopy  (SPECS);  (iii)  galvanos- 
tatic  methods  such  as  chronopotentiometry  and  the  galvanostatic 
intermittent  titration  technique  (GITT);  (iv)  techniques  based  in 
the  frequency  domain  such  as  electrochemical  impedance  spec¬ 
troscopy  (EIS).  To  extract  a  complete  performance  profile  of  the 
material  being  considered,  i.e.,  over  the  entire  discharge  com¬ 
positional  range,  quite  often  it  is  necessary  to  either  combine 
methods  or  establish  an  iterative  sequence  to  achieve  this  goal. 
For  example,  the  techniques  of  LSV/CV  and  chronopotentiom¬ 
etry  are  essentially  used  solely  for  comparative  purposes  since 
the  results  describe  the  kinetic  performance  of  a  given  mate¬ 
rial  only  under  one  set  of  conditions,  i.e.,  voltage  scan  rate  or 
current  density.  The  more  elaborate  SPECS  and  GITT  methods, 
while  providing  both  thermodynamic  and  kinetic  information, 
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are  prone  to  taking  extended  periods  of  time  to  conduct  a  single 
experiment,  in  which  case  material  and  system  stability  becomes 
an  issue.  Similarly,  a  single  EIS  measurement  yields  very  lit¬ 
tle  information  about  the  system  as  a  whole  and  so  it  must  be 
combined  with  another  technique  that  changes  the  solid  elec¬ 
trode  composition  in  a  systematic  fashion.  In  general,  with  solid 
electrode  systems  being  essentially  single-use  only,  i.e.,  past 
electrode  charge-discharge  history  is  important,  it  is  clear  that 
a  trade-off  has  to  be  made  between  the  amount  of  useful  data 
extracted  versus  the  time  taken  to  collect  that  data. 

To  elaborate  further,  SPECS  is  a  very  common  technique  that 
is  known  to  produce  considerable  thermodynamic  and  kinetic 
electrochemical  data  from  a  single  experiment.  A  SPECS  exper¬ 
iment  involves  applying  a  small  voltage  step  (5-25  mV)  to  the 
solid  electrode  that  is  followed  by  an  extended  rest  time  dur¬ 
ing  which  the  resultant  current  is  recorded  as  a  function  of  time 
[1-10].  When  a  pre-defined  experimental  limit  has  been  reached, 

1. e.,  minimal  current  or  a  fixed  time,  the  voltage  is  stepped  again 
and  the  process  repeated  so  that  it  ultimately  covers  the  usable 
discharge  or  charge  compositional  range.  From  the  resultant  cur¬ 
rent  transients,  kinetic  information  such  as  a  diffusion  coefficient 
(as  a  function  of  composition)  can  be  extracted  provided  a  sat¬ 
isfactory  model  is  adopted.  It  is  important  at  this  point  to  be 
clear  that  extraction  of  a  diffusion  coefficient  is  only  appropriate 
when  dealing  with  solid-state  systems  that  undergo  single-phase 
(solid  solution  or  homogeneous)  reduction.  Furthermore,  if  the 
rest  time  is  sufficient,  quasi-thermodynamic  data  can  also  be 
extracted.  While  this  technique  produces  a  large  amount  of  very 
useful  data,  the  time  required  for  the  process  is  a  major  drawback. 
Depending  on  the  conditions  used,  a  SPECS  experiment  can  take 
over  2  weeks  to  complete.  This  not  only  makes  the  technique 
susceptible  to  any  chemical  instability  that  the  electrode  might 
suffer,  but  also  limits  the  number  of  experiments  that  can  be 
conducted  and  exposes  the  overall  process  to  external  problems 
such  as  power  interruptions. 

In  this  study,  an  existing  electrochemical  technique  is  adopted 
to  a  solid  electrode  system  with  the  objective  of  producing 
sufficient  data  to  characterize  material  performance  within  a 
reasonable  timeframe.  Instead  of  applying  a  sequence  of  small 
voltage  steps  as  in  a  SPECS  experiment,  the  approach  has  been 
to  apply  a  much  larger  voltage  step  (0. 1-0.5  V  from  the  open- 
circuit  voltage  (OCV)),  and  then  interpret  the  resultant  data 
with  a  suitable  chronoamperometric  model.  This  experimen¬ 
tal  technique  is  evaluated  using  the  alkaline  manganese  dioxide 
(y-MnCh)  electrode. 

2.  Experimental 

2.1.  Materials 

The  manganese  dioxide  was  provided  by  Delta  EMD  Aus¬ 
tralia,  Pty  Limited  [  1 1  ] .  It  was  prepared  by  anodic  electrolysis  on 
to  a  Ti  substrate  from  a  hot  (~98  °C)  acidic  (H2SO4)  solution  of 
MnSCL,  in  which  the  PhSCLiMn  ratio  was  ~0.4.  When  suffi¬ 
cient  material  had  been  deposited,  the  Ti  anode  was  removed 
from  the  electrolysis  cell  and  the  deposit  was  mechanically 
stripped  from  the  substrate.  The  product  was  subsequently 


milled  to  a  — 105 -pm  powder,  neutralized  and  washed  to  remove 
any  remaining  electrolyte,  and  then  dried  at  110°C.  Materi¬ 
als  produced  in  this  fashion  are  called  electrolytic  manganese 
dioxide  (EMD).  Using  standard  analytical  techniques  it  was 
determined  that  the  EMD  sample  possessed  the  y-crystal  struc¬ 
ture  (X-ray  diffraction  using  Cu  Ka  radiation)  [3],  had  a  OiMn 
ratio  of  1.96  [12],  a  water  content  of  1.72  wt.%  [13-15],  and  a 
BET  surface  area  (No  adsorption  at  77  K)  of  32  m2  g-1  [16], 

2.2.  Electrochemical  cell  construction  [1,2] 

The  working  electrode  blackmix  was  prepared  by  lightly 
grinding  together  in  a  mortar  and  pestle  the  EMD  (0.45  g),  Tim- 
cal  SFG6  graphite  (4.5  g)  and  9  M  (37%)  KOH  (1.105  g).  After 
~5  min  of  mixing,  the  blackmix  was  stored  in  an  airtight  con¬ 
tainer  for  24  h  to  allow  equilibration  to  occur,  i.e.,  the  electrolyte 
to  fully  wet  the  EMD  surface. 

The  electrochemical  cell  is  shown  schematically  in  Fig.  1 .  An 
amount  of  working  electrode  blackmix  corresponding  to  0.015  g 
of  EMD  was  placed  in  a  Teflon-lined  gold-coated  C-size  battery 
can.  The  inner  surface  of  the  Teflon  was  brushed  down  to  remove 
attached  particles  and  three  layers  of  separator  paper  were  laid  on 
top  of  the  sample.  The  mixture  was  then  compressed  under  1 1  for 
60  s  to  form  a  compact  electrode  that  was  in  good  contact  with  the 
can  base.  After  compression,  a  perforated  Perspex  separator  disc 
was  added,  followed  by  a  stainless-steel  counter  electrode.  The 
chamber  was  then  filled  with  ~15  mL  of  9  M  KOH  electrolyte 
and  a  Perspex  cap  was  inserted  into  the  top.  The  cell  was  mounted 
on  top  of  a  brass  current-collector  between  a  Perspex  cover  and 
base  plate,  and  secured  in  place  by  three  bolts  that  were  tightened 
to  a  common  torque  of  0.75  cNm  to  ensure  uniform  pressure.  To 
complete  the  cell,  a  Hg|HgO  reference  electrode  (9M  KOH) 
was  inserted  through  the  Perspex  cap. 

2.3.  Discharge  protocol 

To  determine  the  starting  potential,  the  OCV  of  the  EMD 
was  determined.  Sufficient  EMD  to  form  a  thin  layer  on  the 


Fig.  1 .  Schematic  of  electrochemical  cell. 
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bottom  of  a  plastic  sample  vial  (50-mL  volume)  was  covered 
with  ~20  mL  of  9  M  KOH  solution,  shaken  and  allowed  to  settle 
and  equilibrate  in  a  sealed  environment  for  24  h.  The  OCV  was 
measured  by  placing  a  Pt  disc  electrode  into  contact  with  the 
EMD  in  the  vial  and  measuring  its  potential  with  respect  to  a 
Hg|HgO  reference  electrode  in  the  same  solution. 

Identical  individual  cells  were  prepared  from  the  above- 
mentioned  composite  working  electrode  mixture  for  the 
chronoamperometry  experiments.  This  was  not  intended  to 
reflect  the  proposed  electrochemical  technique,  but  rather 
attempts  to  understand  the  behaviour  of  the  alkaline  EMD 
electrode  under  such  conditions.  Each  cell  was  then  charged 
to  the  OCV  (0.261  V),  at  which  it  was  then  held  for  a  1-h 
rest  period.  The  last-mentioned  step  was  necessary  because, 
in  an  alkaline  environment,  EMD  is  chemically  reduced  by 
the  surface  functional  groups  on  the  graphite.  After  this  rest 
period,  the  potential  was  stepped  cathodically  to  a  predetermined 
voltage  where  it  was  allowed  to  sit  for  7h  with  the  result¬ 
ing  current  recorded  as  a  function  of  time.  The  step  voltages 
employed  were  0.1,  0.0,  —0.1,  —0.2,  —0.3  and  —0.4  V  versus 
Hg|HgO. 

For  the  purposes  of  comparison,  LSV  and  SPECS  experi¬ 
ments  were  carried  out  on  electrodes  constructed  in  the  same 
fashion.  LSV  was  conducted  from  the  OCV  of  the  material 
to  —0.75  V  using  a  scan  rate  of  either  0.02  or  0.05mVs_1. 
The  protocol  used  for  the  SPECS  experiment  involved  sequen¬ 
tial  10 mV  cathodic  steps  with  a  2h  rest  time  between 
steps. 

3.  Results  and  discussion 

3.1.  Chronoamperometry  on  alkaline  manganese  dioxide 
electrode 

The  chronoamperometric  response  of  the  manganese  diox¬ 
ide  electrode  to  various  step  voltages  is  shown  in  Fig.  2(a).  As 
expected,  a  greater  current  resulted  when  the  voltage  was  stepped 
to  more  negative  values.  From  the  i-t  transient,  it  is  quite  straight¬ 
forward  to  convert  this  data  into  chronocoulometric-  and  hence 
compositional-data  as  shown  in  Fig.  2(b).  Again,  as  expected, 
discharge  of  the  electrode  was  greater  when  the  voltage  was 
stepped  to  lower  values. 

Before  developing  a  model  to  extract  quantitative  informa¬ 
tion  from  the  chronoamperometric  data,  considerable  qualitative 
analysis  can  be  carried  out  since  an  i-V-t  relationship  has 
been  established.  To  be  able  to  compare  the  chronoampero¬ 
metric  data  to  linear  sweep  voltammetry,  the  current  at  various 
times  after  the  voltage  step  was  sampled.  By  doing  this  for  all 
of  the  voltage  steps,  a  pseudo-voltammogram  can  be  gener¬ 
ated  from  the  chronoamperometry  data  that  effectively  shows 
voltammetric  discharge  at  different  rates.  Such  data  are  given 
in  Fig.  3,  and  is  compared  with  LSV  experiments  at  two  dif¬ 
ferent  rates.  Interestingly,  there  are  apparent  similarities  in 
the  appearance  of  the  i-V  curves  irrespective  of  how  the  data 
were  collected.  Another  notable  feature  between  the  differ¬ 
ent  data  sets  was  the  time  taken  to  achieve  comparable  data. 
In  Fig.  3(b),  for  instance,  it  appears  that  the  0.02mVs_1 


Fig.  2.  (a)  Chronoamperometric  response  of  alkaline  EMD  electrode;  (b)  cor¬ 
responding  chronocoulometric  and  compositional  data. 


voltammogram  has  a  similar  current  response  to  the  chronoam¬ 
perometric  data  sampled  2000  s  after  the  voltage  step.  This  is 
interesting  given  that  it  took  over  9h  (33  x  103  s)  to  gener¬ 
ate  this  voltammetric  data.  The  similarity  in  current  response 
is  most  likely  explained  simply  by  the  differences  in  the 
techniques  used,  but  it  does  also  suggest  that  equivalent  per¬ 
formance  data  can  be  obtained  in  a  much  shorter  timeframe 
with  the  used  of  multiple  identical  cells  and  chronoamperome¬ 
try. 

In  a  similar  fashion,  the  integrated  form  of  the  chronoam¬ 
perometric  data  (capacity,  Q )  can  be  compared  with  capacity 
data  calculated  from  the  LSV  experiments,  as  shown  in  Fig.  4. 
Again,  the  similarities  between  the  data  sets  are  apparent.  Now, 
however,  the  LSV  data  is  comparable  with  longer  times  after 
the  voltage  step  was  made  in  the  chronoamperometric  experi¬ 
ments,  i.e.,  10-20  x  103  s.  This  is  again  logical  since  the  capacity 
extracted  from  the  electrode  under  these  discharges  at  rela¬ 
tively  low  rates  is  expected  to  at  least  approach  the  theoretical 
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sponding  LSV  data:  (a)  original  view;  (b)  expanded  view. 


capacity  limit  (1109Cg_1),  which  with  chronoamperometry 
experiments  can  only  be  extracted  after  long  periods  of  time 
at  low  voltages. 

The  performance  of  a  battery  electrode  material  can  be  eval¬ 
uated  in  part  by  its  rate  capability  and  capacity,  which  from  a 
strictly  performance  perspective  are  the  two  important  charac¬ 
teristics.  From  the  above  initial  analysis,  it  would  seem  that  large 
voltage  step  chronoamperometry  is  a  valuable  technique  for 
characterizing  both  of  these  electrode  characteristics.  Voltam- 
metric  data  are  related  to  chronoamperometric  data  in  Fig.  3  and 
certainly  from  a  rate  capability  perspective,  the  greater  current 
flowing  at  short  times  after  the  voltage  step,  the  better  is  the  rate 
capability.  With  large  voltage  step  chronoamperometry  the  total 
capacity  is  also  easily  extracted,  here  after  extended  times  at  the 
step  voltage. 


Capacity  /  C/g 

Fig.  4.  Comparison  between  capacity  data  extracted  from  LSV  (0.02  and 
0.05  mV  s-1)  and  chronoamperometry. 


3.2.  Discharge  mechanism  of  alkaline  manganese  dioxide 
electrode 

To  begin  to  model  the  chronoamperometry  data,  it  is  first  nec¬ 
essary  to  elaborate  upon  the  discharge  mechanism  that  the  EMD 
electrode  undergoes  in  an  alkaline  electrolyte.  Many  authors 
have  studied  the  discharge  of  the  alkaline  manganese  dioxide 
electrode  [17-30].  The  basis  of  our  current  understanding  was 
proposed  initially  by  Kozawa  et  al.  [31-35]  who  described  a 
two-step  mechanism.  In  the  first  step  (electron-proton  mecha¬ 
nism),  an  electron  from  the  external  circuit  is  inserted  into  the 
manganese  dioxide  structure  to  reduce  an  Mn4+  ion  to  Mn3+. 
To  maintain  structural  charge  balance,  a  water  molecule  present 
at  the  manganese  dioxide  |  electrolyte  interface  dissociates  into  a 
proton,  which  is  inserted  into  the  structure,  and  an  OH-  ion  that 
remains  in  the  electrolyte.  This  step  in  the  overall  mechanism 
is  a  homogeneous  process  that  leads  to  the  formation  of  a  solid 
solution  of  hydrogen  in  the  host  manganese  dioxide  structure, 


Mn02  +  4H20  +  4e“  -►  Mn00H4  +  ziOH“  (1) 

where  A  is  the  mole  fraction  of  hydrogen  inserted  into  the  struc¬ 
ture  (0<<4<1).  Inserted  protons  diffuse  into  the  bulk  of  the 
structure  by  hopping  to  adjacent  sites  until,  ultimately,  the  hydro¬ 
gen  is  distributed  uniformly  throughout  the  host  structure.  This 
one-electron  reaction  is  primarily  responsible  for  the  power  that 
can  be  extracted  from  the  alkaline  manganese  dioxide  electrode. 
Subsequent  refinements  to  this  step  in  the  mechanism  include 
the  introduction  of  structural  features  in  the  original  manganese 
dioxide  such  as  cation  vacancies  (□),  Mn3+  ions  and  structural 
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water  [13-15],  i.e., 

(Mn4+)1_^_,(n^(Mn3+)y(02-)2-4x-y(OH-)4^ 

+  AH20  +  AsT  -+  (Mn4+)i_x_y_4(D) 

AMn3+)y+A(02-)2-4x-y-A(0H-hx+y+A  +  AOU~  (2) 

where  x  and  y  are  the  mole  fractions  of  cation  vacancies  and 
Mn3+  ions,  respectively.  It  should  be  noted  that  protons  (struc¬ 
tural  water)  are  used  to  compensate  for  charge  discrepancies 
in  the  structure.  A  later  refinement  was  the  identification  of 
structurally  (and  hence  energetically)  different  domains  within 
the  manganese  dioxide  structure,  which  are  hence  reduced  at 
different  voltages  [3,5-10,17-23].  These  different  domains  are 
present  in  the  structure  of  the  manganese  dioxide  produced  by 
electrolysis,  namely  y-MnC>2.  This  polymorph  of  manganese 
dioxide  has  been  described  in  simple  terms  as  a  microscopic 
intergrowth  between  the  pyrolusite  (|3-Mn02)  and  the  rams- 
dellite  forms  of  manganese  dioxide  [36],  as  shown  in  Fig.  5. 
Despite  this  apparent  simplicity,  the  structure  is  complicated  by 
the  cation  vacancies,  Mn3+  ions  and  structural  water  mentioned 
previously,  and  indicated  in  Fig.  5,  as  well  as  by  a  very  finely 
divided  morphology,  which  makes  structural  analysis  difficult. 
Additionally,  micro-twinning  of  the  y-MnC>2  structure  has  been 
postulated  to  account  for  various  features  of  its  X-ray  diffraction 
pattern,  although  the  actual  physical  presence  of  this  feature  has 
still  to  be  confirmed  [3,37], 

The  second  step  in  the  reduction  process,  as  first  described 
by  Kozawa  and  Jeager  [35]  and  refined  by  others  subsequently 
[17-19,38-39],  involves  a  dissolution-precipitation  mecha¬ 
nism.  Here,  the  first-electron  reduction  product  (MnOOH) 
undergoes  dissolution  in  the  concentrated  alkaline  electrolyte  to 
form  a  soluble  Mn(III)  species.  It  has  been  shown  [40,41]  that 
this  mechanism  is  relevant  only  for  more  concentrated  alkaline 
electrolytes  where  the  solubility  of  MnOOH  is  substantial.  Once 
in  solution,  and  at  sufficiently  low  voltages,  the  Mn(III)  species 
undergoes  reduction  at  a  conductive  surface  within  the  cathode 
to  form  a  soluble  Mn(II)  species,  which  precipitates  as  Mn(OH)2 
almost  immediately.  Given  that  Mn(OH)2  is  a  poor  conductor, 
this  causes  passivation  throughout  the  electrode.  Overall,  the 
second  step  in  the  mechanism  occurs  at  too  low  a  voltage  to  be 


(a)  Pyrolusite  (p-MnCh)  (b)  Ramsdellite 


Fig.  5.  Microscopic  intergrowth  of  (a)  pyrolusite  (p-MnCh)  and  (b)  ramsdellite 
from  7-M11O2. 


useful,  and  also  has  poor  efficiency  as  a  result  of  the  resistive 
Mn(OH)2  [17-19]. 

3.3.  Development  of  chronoamperometry  model 


One  of  the  most  common  models  used  to  interpret  chronoam¬ 
perometry  data  is  the  Cottrell  equation  [42],  i.e., 


m  =  m  m 


nFADl'2C* 


(3) 


where  i(t)  is  the  current  flowing  as  a  function  of  time  t  (the 
subscript  d  indicates  diffusion-limited  conditions);  n  is  the  sto¬ 
ichiometric  number  of  electrons  involved  in  the  redox  reaction; 
F  is  the  Faraday  constant;  A  is  the  electrochemically  active  area; 
D  is  the  active  species  diffusion  coefficient;  C*  is  the  bulk  con¬ 
centration  of  active  species.  The  Cottrell  equation  was  derived 
based  on  semi-infinite  planar  diffusion,  and  assuming  that  the 
voltage  step  was  sufficiently  large  so  that  the  surface  concentra¬ 
tion  of  active  species  is  essentially  zero.  From  Eq.  (3),  a  plot  of 
i(t)  versus  f_1/2  should  give  a  straight  line  from  which  the  diffu¬ 
sion  coefficient  can  be  extracted.  Fig.  6(a)  is  such  a  plot  for  the 
chronoamperometry  data  collected  in  this  work.  Clearly  these 
plots  are  not  linear,  for  which  a  number  of  possible  explanations 
may  be  given.  The  first  and  most  obvious  explanation  is  that  the 
boundary  conditions  used  to  derive  Eq.  (3)  do  not  accurately 
represent  the  physical  processes  occurring  in  the  electrode.  This 
is  particularly  true  since  a  solid-state  electrode  is  being  used 
in  which  diffusion  is  an  internal  process  within  each  particle, 
rather  than  into  the  bulk  of  the  electrolyte.  With  a  mean  parti¬ 
cle  size  of  ~45  |xm,  this  means  that  the  semi-infinite  condition 
will  be  unsuitable  not  long  after  the  voltage  step  because  of 
intersecting  diffusion  fronts  from  opposite  sides  of  the  particle. 
Furthermore,  the  assumption  that  planar  diffusion  is  apparent  is 
also  unreasonable  given  the  morphology  of  the  particles,  which 
is  very  irregular  on  both  the  macroscopic  and  microscopic  scale. 
Both  of  these  features  are  also  complicated  quite  considerably 
by  material  porosity,  which  increases  the  active  surface  area  by 
at  least  two  orders  of  magnitude  compared  with  the  geometric 
area  based  on  particle  size.  Finally,  the  crystal  structure  of  the 
active  material  also  plays  a  part  in  determining  the  chronoamper- 
ometric  response.  As  described  above,  and  shown  in  Fig.  5,  the 
structure  of  y-MnCh  contains  energetically  different  domains 
which,  for  the  purposes  of  these  experiments,  means  that  they 
are  reduced  at  different  rates,  i.e.,  different  diffusion  coefficients. 
In  summary,  there  may  be  overlapping  processes  apparent  in 
the  chronoamperometric  data,  the  sum  total  of  which  causes 
deviation  from  the  behaviour  predicted  by  Eq.  (3). 

To  demonstrate  this  using  the  Cottrell  equation,  we  can  calcu¬ 
late  for  each  of  the  voltage  steps  the  value  it112,  which  essentially 
represents  the  collection  of  constants  for  a  particular  reduction 
process.  Such  data  is  shown  in  Fig.  6(b)  in  comparison  with  elec¬ 
trode  composition.  Clearly  it112  (for  the  most  part)  decreases  over 
the  entire  composition  range.  It  is  concluded  that  this  indicates, 
at  least  within  the  bounds  of  the  present  simple  model,  that  either 
(or  both)  A  and  D  decreases  with  composition,  since  C*  should 
increase  with  an  increasing  step  voltage.  The  only  variation  is 
for  the  step  to  the  lowest  voltage  (— 0.4  V)  where  the  influence 
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x  in  MnOx 


Fig.  7.  Comparison  between  experimental  and  modelled  data  based  on  spherical 
diffusion  into  a  solid  electrode. 

adapted  here  to  consider  large  voltage  steps.  This  model  is  based 
on  spherical  diffusion  into  a  solid  electrode  material  using  the 
following  initial  and  boundary  conditions: 

Initial  :  C(x,  t )  =  Co,  for  all  x  when  t=  0  (4) 

Boundary  :  C(x,  t)  =C\,  forx=  0  when  t  >  0  (5) 

B  dC(x,  t )  i(t) 

Boundary  :  - =  - ,  for  x  =  0  and  a 

ox  nFAD 

when  t  >  0  (6) 

Solution  of  Fick’s  second  law  of  diffusion  under  such  condi¬ 
tions  leads  to: 

2FADAC ^  n2n2Dt 

i  =  - V - 2 —  (7) 


Fig.  6.  (a)  Chronoamperometry  data  analysis  based  on  Cottrell  equation  (j-r  1/2) 
and  (b)  it112  as  a  function  of  composition. 

of  the  second  electron  reduction  is  apparent.  A  decrease  in  A  can 
be  explained  by  the  structural  expansion  that  occurs  upon  reduc¬ 
tion  [43] .  With  the  y-MnC>2  being  microporous,  any  swelhng  due 
to  the  inclusion  of  the  larger  Mn3+  ions  into  the  structure  may 
cause  the  micropores  to  close  down  and  thereby  limit  the  active 
surface  area  of  the  electrode  material.  Similarly,  the  drop  in  D 
can  be  justified  by  considering  the  proton  transport  mechanism 
associated  with  y-MnOo,  since  this  is  believed  to  be  the  rate- 
limiting  phenomena  [20].  Proton  transport  occurs  via  hopping, 
and  is  therefore  dependent  on  the  availability  of  unoccupied  sites 
in  the  direction  of  diffusion.  With  a  more  reduced  material,  the 
number  of  unoccupied  sites  decreases,  essentially  inhibiting  the 
rate  of  proton  transport. 

A  much  more  advanced  diffusion  model  was  developed  by 
Hong  et  al.  [22],  albeit  for  small  potential  step  experiments,  but 


where  most  symbols  have  already  been  defined,  except  for  AC 
which  is  the  change  in  concentration  of  hydrogen  in  the  man¬ 
ganese  oxide  structure  over  the  potential  step  (C\  —  Co),  and  a 
is  the  radius  of  the  particle.  By  combining  constants  in  Eq.  (7) 
the  fitting  parameters  P\  and  P2  can  be  produced,  i.e., 

i  =  Pi^exp(-n2P20  (8) 


where 

2FADAC 
Pi  = - 


and 


(9) 

(10) 


Pi  and  P2  can  now  be  determined  using  linear  least-squares 
regression  by  fitting  curve(s)  defined  by  Eq.  (7)  to  the  experimen- 
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tal  data.  The  comparison  between  the  predicted  and  experimental 
chronoamperometric  data  is  shown  in  Fig.  7,  where  an  attempted 
is  also  made  to  address  the  need  for  multiple  curves  defined  by 
Eqs.  (7)  and  (8)  to  account  for  the  reduction  of  energetically 
different  sites  within  the  y-MnCh  electrode.  The  experimental 
data  used  in  this  figure  are  for  the  step  to  —0.3  V.  Fitting  two 
curves  leads  to  a  result  that  quite  closely  matches  the  exper¬ 
imental  data  (average  deviation  of  ±0.004  A  g-1).  A  similar 
deviation  is  apparent  for  all  other  data  sets  that  represent  steps 
to  other  voltages.  Fitting  a  greater  number  of  curves  (>2)  to 
the  experimental  data  does  not  improve  the  quality  of  the  fit 
significantly  (statistically,  F-test)  enough  to  warrant  their  inclu¬ 
sion,  despite  the  already  demonstrated  presence  of  more  than 
two  processes  [1,2],  For  two  fitted  curves,  the  greatest  devia¬ 
tion  between  the  fitted  and  experimental  data  occurs  at  times 
less  than  ~2  x  103  s,  which  represents  the  transition  from  high- 
rate  discharge  to  a  more  sustained  low-rate  extraction  of  charge. 
While  this  definition  is  somewhat  arbitrary,  it  does  perhaps  repre¬ 
sent  complete  utilization  of  certain  domains  within  the  structure, 
while  others  are  continuing  to  be  reduced  at  a  much  slower  rate. 
This  is  also  reflected  in  the  inset  data  in  Fig.  7  where  there  is 
a  very  fast  process  in  parallel  with  a  much  slower  discharge 
process.  It  is  also  apparent  that  this  transition  is  the  greatest 
source  of  error  in  the  modelling.  P\  and  P2  for  each  of  the 
two  processes  are  shown  in  Table  1  as  a  function  of  the  step 
voltage. 

P\  and  Pi  can  also  be  used  to  calculate  A y/D  for  each  process, 
as  shown  in  Eq.  (11): 


A-/D  = 


[Pf 

2FACy  P2 


(11) 


This  expression  is  derived  by  eliminating  the  particle  radius  (a) 
from  Eqs.  (9)  and  (10),  such  that  its  value  does  not  need  to  be 
known.  The  significance  of  Ay/D  is  that  it  is  an  ideal  summation 
for  the  kinetics  of  an  individual  discharge  process,  since  it  avoids 
the  need  to  assume  a  value  for  the  electrochemically  active  sur¬ 
face  area,  which  for  porous  electrode  materials  is  not  known 
with  a  great  degree  of  certainty.  As  mentioned  above,  AC  is  the 
change  in  hydrogen  concentration  in  the  manganese  structure 
during  discharge.  For  each  process  apparent  in  the  chronoamper¬ 
ometric  discharge,  the  first  step  in  determining  AC  is  to  calculate 
the  total  charge  passed  as  a  result  of  the  voltage  step.  From 
the  X-ray  diffraction  pattern  for  the  material  (Fig.  8),  it  is  then 
possible  to  determine  the  unit-cell  dimensions  of  the  y-Mn02 


I 


E 

I 
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Fig.  8.  X-ray  diffraction  pattern  of  EMD  together  with  peak  assignments  based 
on  orthorhombic  unit  cell. 

assuming  orthorhombic  symmetry  [3]  and  four  formula  units 
per  unit  cell  [44],  The  charge  passed  and  unit-cell  dimensions 
then  allow  determination  of  the  average  number  of  hydrogen 
atoms  per  unit  volume  (concentration),  and  hence  the  change 
that  occurs  with  the  voltage  step.  The  assumptions  made  with 
this  calculation  are  that  y-MnCh  can  indeed  be  described  by  an 
orthorhombic  unit  cell  (given  the  tetragonal  nature  of  (3-Mn02 
which  is  a  structural  component),  and  that  the  unit-cell  volume 
remains  constant  with  depth-of-discharge.  It  is  recognized  that 
this  last  assumption  is  not  strictly  true  since  there  have  been 
a  number  of  previous  reports  that  indicate  the  structure  swells 
on  discharge  due  to  the  inclusion  of  protons  and  Mn(III)  in  the 
structure  [43].  The  expansion  is  at  most  18%  based  on  a  com¬ 
parison  of  the  unit-cell  dimensions  of  ramsdellite  (MnCh)  and 
groutite  (a-MnOOH)  as  reported  in  the  literature  [44] .  While  this 
expansion  is  significant,  the  actual  relationship  between  unit-cell 
volume  and  depth-of-discharge  is  not  linear,  given  the  structural 
complexities  of  the  material.  Hence,  this  is  the  best  estimate  at 
this  time. 

An  important  point  to  mention  at  this  stage  is  that  Ay/D 
cannot  be  separated  into  individual  A  and  D  values  because 
the  electrochemically  active  surface  area  for  manganese  diox¬ 
ide  is  not  known  with  certainty.  The  reason  for  this  is  due  to 


Table  1 

Data  extracted  from  modelling  chronoamperometry  data 


Step  voltage 

Process  1 

Process  2 

Pi 

Pi 

Qt  (C  g-1) 

Ay/D  (x  1010  m3 

s  1/2  g  ') 

Pi 

Pi 

Qi  (Cg-1) 

Ay/D(x  1010  m3  s-1/2  g-1) 

0.1 

0.141 

6577 

18 

4.92 

0.022 

469 

45 

1.17 

0.0 

0.402 

3897 

92 

3.64 

0.011 

59 

170 

0.42 

-0.1 

0.668 

6212 

93 

4.76 

0.050 

89 

544 

0.51 

-0.2 

0.826 

3720 

198 

3.56 

0.107 

161 

644 

0.68 

-0.3 

1.607 

5639 

247 

4.48 

0.237 

326 

702 

0.97 

-0.4 

0.880 

1276 

645 

2.00 

0.035 

74 

460 

0.46 
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the  discrepancy  that  exists  between  the  methods  that  are  used 
to  measure  surface  area  and  the  species  that  are  electrochemi- 
cally  active  when  manganese  dioxide  is  immersed  in  an  alkaline 
electrolyte.  The  most  common  method  for  surface  area  mea¬ 
surement  is  N2  gas  adsorption  at  77  K,  which  is  subsequently 
interpreted  using  the  BET  isotherm  [16].  Typical  BET  surface 
areas  for  commercial  EMD  lie  in  the  range  of  28-32  m2  g_1, 
although  EMD  samples  have  been  produced  with  BET  surface 
areas  that  range  from  10  to  100  m2  g-1.  Overall,  these  values 
indicate  that  EMD  is  a  very  porous  material  and  that  in  actual 
fact  the  surface  area  is  mostly  associated  with  micropores  [45]. 
Compare  this  method  of  analysis  with  the  alkaline  manganese 
dioxide  electrode  interface  in  which  the  H2O  molecule  is  the 
important  adsorbate.  Clearly  there  are  significant  differences  in 
molecular  size,  shape,  polarity,  etc.,  that  really  limit  the  appli¬ 
cability  of  N2  adsorption  as  a  characterization  method.  Despite 
the  apparent  differences  in  adsorbates,  we  also  need  to  differen¬ 
tiate  between  surface  area  and  pore  size.  As  mentioned  above, 
EMD  is  a  microporous  material,  and  on  this  scale  understand¬ 
ing  the  mass  transport  of  water  into  and  out  of  the  EMD  pores  is 
also  critical.  Because  water  is  consumed  in  the  electrochemical 
reaction  to  form  OH-  ions,  water  must  diffuse  from  the  bulk 
solution  into  the  pore  to  replenish  the  amount  that  is  consumed, 
and  at  the  same  time  OH-  ions  must  diffuse  out  of  the  pore  to 
prevent  surface  passivation.  Therefore,  the  rate  of  electrochem¬ 
ical  discharge  will  determine  the  active  surface  area  since  mass 
transport  of  electroactive  species  within  material  pores  exerts  an 
effect. 

The  calculated  values  of  Ay/D  for  each  of  the  two  pro¬ 
cesses  are  also  shown  in  Table  1  as  a  function  of  step  voltage. 
Clearly,  Process  1  is  kinetically  faster  (~4  times)  than  Process 
2.  This  calculation,  coupled  with  the  much  smaller  capacity 
associated  with  Process  1  (for  the  most  part),  indicates  that 
it  represents  a  structural  domain  that  is  readily  consumed  at 
a  fast  rate.  The  origin  of  this  domain  is  not  known  with  cer¬ 
tainty,  but  it  may  be  related  to  reduction  of  Mn(IV)  ions  near 
defects  such  as  cation  vacancies  in  the  structure.  Furthermore, 
our  previous  work  [2]  has  shown  that  reduction  of  pyrolusite 
domains  within  the  structure  is  also  a  kinetically  fast  process. 
With  regards  to  material  design,  this  represents  a  desired  domain 
within  the  y-MnCb  structure  for  high-rate  performance  materi¬ 
als.  Process  2,  on  the  other  hand,  may  represent  the  slightly 
slower  reduction  of  the  majority  ramsdellite  domains  within  the 
structure.  Overall,  the  values  obtained  here  are  consistent  with 
those  reported  previously  for  studies  on  y-Mn02 .  Many  attempts 
have  been  documented  in  the  literature  concerning  the  determi¬ 
nation  of  proton  diffusion  coefficient  in  y-MnCE.  A  wide  range 
of  values  have  been  reported  (10-11-10-18  cm2  s-1)  based  on 
measurements  made  using  a  variety  of  methods  [46].  If  these 
measurements  of  D  are  combined  with  the  area  assumed  in  the 
calculation,  then  Ay/D  irrespective  of  the  method  can  be  calcu¬ 
lated  to  be  within  an  order  of  magnitude  of  ~  10“  3  cm3  s-1/2  g~  1 . 
The  work  reported  here  is  consistent  with  this;  it  gives  an 
average  value  of  3.89  x  10_4cm3  s-1/2g-1  for  Process  1  and 
0.70  x  10-4  cm3  s_1/2  g-1  for  Process  2.  Aside  from  taking  an 
average  value  for  Ay/D,  there  does  not  appear  to  be  any  sig¬ 
nificant  trends  in  the  data  with  step  voltage,  which  is  expected 


step  chronoamperometry. 

given  that  the  same  processes  are  being  examined  throughout 
the  discharge  irrespective  of  the  step  voltage. 

The  value  of  the  large  voltage  step  chronoamperometry 
experiments  becomes  apparent  when  the  resultant  data  are  com¬ 
pared  with  a  SPECS  experiment  on  the  same  electrode.  In  this 
case,  a  10-mV  cathode  step  was  applied  every  2  h  from  the  open- 
circuit  voltage  to  less  than  —0.4  V.  A  similar  approach  was  used 
to  model  the  i-t  transients,  i.e.,  the  spherical  model  described 
by  Eq.  (7)  combined  with  the  subsequent  analysis  to  determine 
Ay/D  (Eq.  (11)),  taking  into  proportionate  account  the  multi¬ 
ple  processes  occurring  at  different  voltages  (as  described  in 
Ref.  [2]),  and  the  fact  that  only  a  small  cathodic  voltage  step 
was  used.  The  resultant  A^/D  data  is  shown  in  Fig.  9  as  a 
function  of  electrode  voltage  (and  hence  composition),  where 
it  is  compared  with  the  results  obtained  using  large  voltage  step 
chronoamperometry.  For  the  SPECS  data,  peaks  A—F  refer  to  the 
component  processes  that  were  extracted  in  the  SPECS  exper¬ 
iment  from  the  overall  y-Mn02  discharge  process.  There  is  a 
good  comparison  between  the  Ay/D  values  obtained  using  the 
different  methods,  thus  establishing  the  validity  and  applicabil¬ 
ity  of  our  large  voltage  step  chronoamperometry  method.  As 
mentioned  before,  Process  1  reflects  the  higher  discharge  rate 
processes,  while  Process  2  represents  the  more  sustain,  longer 
lasting  extraction  of  charge. 

4.  Summary  and  conclusions 

The  applicability  of  large  voltage  step  chronoamperometry  to 
the  study  of  solid  electrodes  has  been  evaluated,  with  the  alkaline 
y-MnCK  electrode  representing  a  test  case.  The  results  obtained 
are  compared  with  those  of  more  conventional  LSV  and  SPECS 
analysis  methods.  The  motivation  of  the  work  is  to  develop 
an  alternative  electrochemical  characterization  method  that  is 
quicker  but  still  provides  valuable  information.  While  steps  to 
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multiple  voltages  are  time-consuming  experiments  when  com¬ 
bined  they  provide  the  ability  to  evaluate  the  method.  Routine 
analysis  using  this  method  is  to  only  one  step,  voltage,  namely 
—0.3  V,  so  as  to  examine  almost  all  of  the  one-electron  discharge 
and  avoid  complications  associated  with  accessing  reduction 
with  the  second  electron. 

From  an  initial  qualitative  analysis,  it  has  proved  possible  to 
demonstrate  similarities  between  the  chronoamperometry  and 
LSV  performance  data.  In  particular,  by  sampling  the  current 
at  defined  times  after  the  voltage  step,  data  have  been  gener¬ 
ated  that  are  comparable  with  the  raw  and  integrated  LSV  data 
(Figs.  3  and  4,  respectively). 

The  chronoamperometry  data  are  then  modelled  using 
Fick’s  second  law  and  assuming  spherical  diffusion  into  solid 
particles.  Multiple  curves  (2)  are  necessary  to  match  the 
experimental  data,  with  two  curves  that  statistically  represent 
fast  and  slow  processes.  From  this  analysis  values  for  A^/D 
are  extracted  (3.89  x  10_4cm3  s_1/2g_1  for  Process  1  and 
0.70  x  10_4cm3  s-1/2g-1  for  Process  2)  are  found  to  match 
quite  closely  those  reported  previously  determined  by  SPECS 
using  the  much  more  time-consuming  and  complex  SPECS  tech¬ 
nique  (Fig.  9). 

In  conclusion,  the  effectiveness  of  large  voltage  step 
chronoamperometry  as  an  efficient  method  for  characterizing 
solid  electrode  materials  has  been  demonstrated. 
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Glossary  of  terms 

General  terms 

CV:  Cyclic  voltammetry 
EIS:  Electrochemical  impedance  spectroscopy 
GITT:  Galvanostatic  intermittent  titration  technique 
LSV:  Linear  sweep  voltammetry 
OCV:  Open-circuit  voltage  (V  or  mV) 

Processes  1  and  2:  Individual  discharge  processes  during  the  overall  electrode 
discharge 

SPECS:  Step  potential  electrochemical  spectroscopy 
Chronoamperometric  model  terms 
a:  Spherical  particle  radius  (m  or  cm) 

A:  Electrochemically  active  surface  area  (m2  g-1  or  cm2  g-1) 
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C*:  Bulk  concentration  (mol  L-1  ormolm-3)  i(t):  Time-dependent  current  (A  g-1  ormAg-1) 

Co:  Initial  concentration  (mol  L-1  ormolm-3)  i<i(t):  Diffusion-limited  current  (A  g-1  ormAg-1) 

C(x,  t):  Spatial  and  time  dependent  species  concentration  n:  Stoichiometric  number  of  electrons  transferred 

A  C:  Change  in  species  concentration  after  discharge  (mol  L- 1  or  mol  m-3 )  Q:  Capacity  (C  g- 1 ) 

D:  Diffusion  coefficient  (m2  s-1  or  cm2  s-1)  t:  Time  (s) 

F:  Faraday  constant  (96,485  C  mol-1) 


